Semiclassical theory of the quasi two-dimensional trapped Bose gas 
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We discuss the quasi two-dimensional trapped Bose gas where the thermal occupation of excited 
states in the tightly confined direction is small but remains finite in the thermodynamic limit. 
We show that the semiclassical theory describes very accurately the density profile obtained by 
Quantum Monte Carlo calculations in the normal phase above the Kosterlitz-Thouless temperature 
Tkt, but differs strongly from the predictions of strictly two-dimensional mean-field theory, even at 
relatively high temperature. We discuss the relevance of our findings for analyzing ultra-cold-atom 
experiments in quasi two-dimensional traps. 
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For many years, the physics of two-dimensional quan- 
tum gases has been under close experimental and theoret- 
ical scrutiny. The quest for quantum phase transitions in 
two-dimensional atomic Bose gases has started with ex- 
periments on spin-polarized atomic hydrogen adsorbed 
on liquid '*He surface reaching the quantum degenerate 
regime [l| and observing the onset of quasi-condensation 
0. The Kosterlitz-Thouless transition Q was observed 
recently in trapped atomic gases of ultra-cold ^^Rb atoms 
in an optical lattice potential with a tightly confined z- 
direction 

In contrast to experiments with liquid ^He films 
where the Kosterlitz-Thouless transition is realized di- 
rectly, the typical extension in the ^-direction in the ex- 
periments on two-dimensional gases is much larger than 
the three-dimensional scattering length. For this rea- 
son, the effective two-dimensional interaction strength re- 
mains sensitive to the density distribution in z Q . Never- 
theless, the gas is kinematically two-dimensional because 
of the strong out-of-plane confinement. 

In this letter, we consider the quasi two-dimensional 
regime of the trapped Bose gas, where the temperature 
T is of the order of the level spacing in z. This corre- 
sponds to the experimental situation with small, but not 
completely negligible, thermal occupation of a few ex- 
cited states of the tightly confining potential 0, Q. The 
quasi two-dimensional regime crosses over to the three- 
dimensional and the strictly two-dimensional Bose gases 
as the potential in the z-direction is varied. 

It was noticed in experiments^, Q and in direct Quan- 
tum Monte Carlo calculations [9| that the density profile 
of the gas deviates strongly from two-dimensional mean- 
field theory, even at relatively high temperature. We 
point out that these deviations can be accounted for by 
a quasi two-dimensional mean-field theory which incor- 
porates corrections due to the tightly confined direction. 

We first discuss Bose-Einstein condensation of the 
ideal gas in strongly anisotropic harmonic traps. We de- 



fine the quasi two-dimensional regime where the ideal- 
gas critical temperature is always lower than that of the 
strictly two-dimensional ideal gas. Including interactions 
on the level of mean-field, we obtain the density pro- 
files in the semiclassical approximation and solve the 
self-consistent mean-field equations directly. Remark- 
able agreement of the semiclassical density profiles with 
the results of Quantum Monte Carlo calculations is ob- 
tained in the high-temperature normal phase down to the 
Kosterlitz Thouless temperature. The profiles should be 
very convenient for calibrating the temperature in exper- 
iments of quasi two-dimensional Bose gases. Compari- 
son of experimental density profiles with Quantum Monte 
Carlo data has already removed the original discrepancy 
of the Kosterlitz-Thouless temperature between calcula- 
tion and experiment 0]. 

We consider an anisotropic trap with oscillator fre- 



<^ ujz- At temperature T 



y and 



quencies to = uj^ = 

the motion is semiclassical in the coordinates 
in the momenta Hk^, Hky, whereas the quantization in 
the 2-direction is best described through the energy lev- 
els lyHiOz {j^ = 0, 1, ... ) of the corresponding harmonic 
oscillator. Semiclassically, the number AN of particles 
per phase space element dkxdkydxdy in the energy level 
V is given by [lol | 



dN = 



dkxdkydxAy 



(27r)2 cxp + v{r) + vhio, - - I 



, (1) 



where (5 = 1/7", k"^ = k'^ + ky, and where v{r) is an 
arbitrary two-dimensional potential energy (with = 
x^ + y^). 

Equation ([!]) can be integrated over all momenta and 
summed over all oscillator levels to obtain the two- 



2 



dimensional particle density 



n{r) = ^ / 



= - ^ £ ln{l - exp [/3(/x - v{r) - i^hw,)]}, (2) 



where A = y/2TTh?/3/m is the thermal wavelength. The 
potential v{r) can itself contain the interaction with the 
density n{r), so that eq. ^ is in general a self-consistency 
equation. The integral of n{r) over space yields the equa- 
tion of state, that is, the total number of particles as a 
function of temperature and chemical potential. 

Let us first consider the ideal gas, where the potential 
energy v{r) = is due only to the trapping po- 

tential, so that the rhs of eq. ([2]) is independent of the 
density n{r). We get 



N 



i/=0 



1 - e 



0{fi — i'huj~ —muP'r^ 12') 



ft2 



(3) 



where we have defined 



n=l 



The saturation number TVgat (T) is the maximum num- 
ber of excited particles (reached at ^ = 0) at a given 
temperature. We have 



rp2 



1^=0 



The above relation between the saturation number and 
the temperature defines the dependence of the Bose- 
Einstein condensation temperature on the particle num- 
ber N. As mentioned before, the strictly two-dimensional 
limit is characterized by the limit phiOz — > oo (the level 
spacing in z is much larger than the temperature). In 
this limit, only the first term in eq. ([3]) contributes. Us- 
ing ^2(0) = 7r^/6, we find 



h?LL)'^ 6 



(5) 



In the quasi two-dimensional case, with finite Phw^, the 
occupation of the oscillator levels v = 1,2,... increases 
the saturation number and therefore lowers the critical 
temperature. It is convenient to express in units of Tb^q 
both the temperature t = T/T^'^q and the oscillator 
strength Cjz = fvjJz/T^Q, and to write jl = /3/i. Us- 
ing eq. ([S|), we may rewrite the equation of state, eq. ([3]), 




FIG. 1: Bose-Einstein condensation temperature fBEC of 
the ideal quasi two-dimensional gas in a harmonic trap with 
cjz/oj oc A''^''^ (expansions from eq. (|8])). The inset shows tBEC 
for scaling uj^/uj = A^** as a function of 5 for different A''. At 
5 = 0, the critical temperature TbIc = [^/C(3)]^'^^ of the 
ideal gas in a three-dimensional isotropic trap is recovered. 



as a relation between the temperature i, the chemical 
potential /i, and the oscillator strength 



with 



i'=0 / 



-1/2 



(6) 



(7) 



Equation ([6]) is solved numerically by iterating tn+i = 
f{tn) from an arbitrary starting temperature to to the 



(4) fixed point. The critical temperature ^bec = Tbec/-^] 



n2d 
^BEC 

(as a function of Cjz) of the quasi two-dimensional ideal 
Bose gas is the solution for /i = (see fig. [T|) . The re- 
duction with respect to the strictly two-dimensional case 
is notable for systems of experimental interest. For ex- 
ample, we find iBEC = 0-78 for the experimental value 
Cbz = 0.55 [3,[3| considered in the Quantum Monte Carlo 
calculations [9|. 

We can expand eq. d?]) for small and for large Uz and 
find 



*BEC 



02) 
C(3) 



1/3 



-1/3 



1C(2)~, 
6 C(3) 



for (D, ^ 1 



2C(2)-^ 



■ exp {—ujz) for ujz^ 1 



(8) 



where we have used C(s) = -Ps(O) (note that C,{2) = 77^/6 
and C(3) ~ 1.202). The expansions are indicated in 
fig. [1] They give the critical temperature to better than 
1.2% for all values of uJz (the low-ujz expansion is used 
for uj < 1.8 and the high-diz expansion for lu > 1.8). 
The first term in the small-w^ expansion of eq. ^ cor- 
responds to the three-dimensional gas. Indeed, iBEC ~ 



3 



[C(2)/C(3)]^/^<I'y^ is equivalent to 



us to determine the equation of state explicitly: 



Tbec 



C(2) 



1/3 



C(3) 



(9) 



the well-known condensation temperature for the three- 
dimensional Bose-Einstein gas in an anisotropic trap [11 1 . 
Equation ^ also follows directly from eq. ^ by replac- 
ing the sum over the oscillator levels by an integral. Like- 
wise, the first term in the large-w^ expansion of eq. ([8]) 
represents the strictly two-dimensional gas. 

The inset of fig. [1] further analyzes the expansions of 
eq. ([8]). Indeed, we can choose a scaling LOzjio ^ dif- 
ferent from the quasi two-dimensional case S = 1/2. Any 
choice oi S < 1/2 corresponds to (D^ — > for iV ^ oo 
so that asymptotically the three-dimensional regime is 
reached. Analogously, i5 > 1/2 corresponds to ^ oo 
for iV — > oo, driving the transition into the strictly two- 
dimensional regime. The rescaled transition tempera- 
tures are plotted for lUz/uj = where the case 5 = 
corresponds to the three-dimensional isotropic trap [l2| • 

In order to describe interaction effects in the quasi two- 
dimensional Bose gas, we now add a semiclassical contact 
term to the potential energy of eq. ([2]): 



v{r) = + 2g[n{r) - ?i(0)]. 



(10) 



We have subtracted the central density, so that the po- 
tential still vanishes at the origin. As discussed earlier 
0, the effective interaction g = ATrafi^/m J dz [p{z)]'^ is 
proportional to the three-dimensional s-wave scattering 
length a and to the integral of the squared density distri- 
bution in z, described by the normalized diagonal den- 
sity matrix p{z). In the temperature range of interest, 
this density distribution is well described by the single- 
particle harmonic-oscillator density matrix in z, leading 
to 



g = a\\ Vtanh[cc)2/(2<)J. 



(11) 



The effective interaction thus decreases with tempera- 
ture from its zero-temperature value g ~ a^^-Kuj^fi^ jm. 
To keep the interaction strength fixed, we must keep g 
(or equivalently g) constant in the quasi two-dimensional 
thermodynamic limit which requires a fixed value of the 
scaled scattering length a^uJl oc au^/'^N^/^. 

The mean- field density 7i(r), on the Ihs of eq. ([2]), de- 
pends on the variable r only via the potential v{r). In 
the space integral over the particle density, we can thus 
change the integration variable from r to v. This allows 
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ov 
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(12) 



where the central density, 



n(0)A2 



i/=0 



In { 1 — exp [fl 



(13) 



is directly expressed in terms of /2, independent of the 
interaction, due to the subtraction performed in our ef- 
fective potential, eq. ((T0|) . (Note that the first integral on 
the second line of eq. ([T^ has already appeared in eq. ^ 
and that the second integral is a total derivative.) The 
equation of state can be written in terms of the tempera- 
ture t. This yields the following generalization of eq. ([7]) 
to the mean-field gas: 
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mg 
2^ 



[n{0)X 



2l2 



-1/2 



(14) 



An iteration procedure 

■^BEC 



f{tn) again obtains the 
temperature t — T/T^q as a function of the chemical 
potential jl for given parameters Coz and g. (The central 
density is computed using eq. (jl3p during each iteration.) 

The mean-field density profile is obtained in two steps 
by first calculating the density profile as a function of the 
scaled effective potential v = (3v, 

oo 

n{v))? = — In {1 — exp (/i — u — lyuiz/t)} , 

i/=0 

and then by inverting eq. (llOp to obtain r(n) (thus n(r)) 
for the given v and nX^: 



r{n}? ,v) 



2T 



mgnX^ 



(15) 



In fig. [21 we show the remarkable agreement of the 
quasi two-dimensional mean-field profile with the one 
obtained by Quantum Monte Carlo simulations as in 
Ref. 0, for N = 10 000 bosons for parameters ujz ~ 0.55 
t = r/T||c = 1, and mg/fi^ = 0.13. The simula- 
tions take into account the full three-dimensional geom- 
etry, and particles interact via the three-dimensional s- 
wave scattering length a (see Ref. [l^l for a more de- 
tailed description of finite temperature simulations of 
trapped Bose gases). Comparison with the ideal quasi 
two-dimensional gas is also very favorable. 
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FIG. 2: Two-dimensional density profile n{r = i/aj^ + y^)\^ 
at temperature T = Tbec in ^ trap with = O.SSTbeci in 
the ideal Bose gas and for rrig/h^ = 0.13 according to mean- 
field theory, compared to Quantum Monte Carlo simulations 
at A'^ = 10 000, and for ideal distinguishable particles (first 
term in eq. Uni). The inset shows the cycle weights nk for 
the strictly two-dimensional and the quasi two-dimensional 
Bose gas, and for the Quantum Monte Carlo simulations (from 
above). 



FIG. 3; Two-dimensional density profile n{r — i/aj^ -f y^)X^ 
at temperature T = 0.8 Tbec in a trap with lu^ = O.SSTbeci 
in the ideal Bose gas and for mg/h^ = 0.13 according to 
strictly two-dimensional and quasi two-dimensional mean- 
field theory (using g and g), compared to Quantum Monte 
Carlo simulations with A*" = 10 000. The inset shows the 
Kosterlitz-Thouless temperature within modified mean-field 
theory as a function of Cj^ ( the strictly two-dimensional limit 
is tKT = 0745 (from eq. 



The mean-fiGld density is an exact sum of Gaussians, 



N 



n(r)A^ 



^ ^ klTk cxp 



A;=l 



(16) 



whose variances correspond to the density distribution 
of the harmonic oscillator at temperature The pre- 
factors in cq. (fT6|) contain the cycle weights TTfc. These 
weights give the probability of a particle to be in a cy- 
cle of length k in the path-integral representation of the 
Bose gas, where the density matrix must be symmetrized 
through a sum of permutations 3, 13] ■ For the ideal Bose 
gas, the TTfe are easily computed. In the distinguishable- 
particle limit, at infinite temperature, only cycles of 
length fc = 1 contribute (tti = 1), whereas the Bose- 
Einstein condensate is characterized through contribu- 
tions of cycles of length k oc N. The cycle weights tt^ are 
shown in the inset of fig. [2] for small k. Only very short 
cycles contribute, and the density profile can thus be de- 
scribed by a very small number of Gaussians. The exact 
cycle-weight distribution of the Quantum Monte Carlo 
does not rigorously correspond to a profile as in cq. (fT6|) . 
although the corrections are negligible in our case. 

Figure [3] considers the temperature t = 0.8 in the inter- 
val between the Kosterlitz-Thouless temperature of the 
quasi two-dimensional interacting gas (at t = txx — 0.70 
for these parameters and the strictly two-dimensional 
Bose-Einstein condensation temperature. Again, the 
agreement of the quasi two-dimensional mean field with 
the exact density profile obtained by Quantum Monte 
Carlo is remarkable. At this temperature, the deviations 



with the ideal quasi two-dimensional profile and with the 
strictly two-dimensional mean field are important. To 
illustrate the temperature-dependence of the effective in- 
teraction, we also show the mean-field profile computed 
with the zero-temperature interaction parameter g in- 
stead of the true effective two-dimensional interaction g 
(see eq. (fTTjl ). We note in this context that the differ- 
ence in eq. (|lip between g and g was determined under 
the condition that the interaction leaves the density dis- 
tribution in z unchanged. Whenever this condition is 
violated, the density distribution in z must be computed 
by other means, as for example by Quantum Monte Carlo 
methods (see [9[). 

Let us finally discuss the Kosterlitz-Thouless transi- 
tion into the low-temperature phase, which is not con- 
tained in mean-field theory. The semiclassical quasi two- 
dimensional gas does not Bose-condense because the par- 
ticle number in eq. (|12|) diverges at p, = (it saturates at 
a finite value in the ideal Bose gas). This divergence is 
due to the logarithmic divergence of the central density 
(see eq. (|13p ). However, interaction effects beyond mean- 
field drive a Kosterlitz-Thouless phase transition Q from 
the high-temperature normal phase to a superfiuid below 
Tkt- 

As discussed previously 0, [3], the Kosterhtz- 
Thouless transition occurs when the central density 
ri(0)A^ reaches the critical value of the two-dimensional 
homogeneous gas, which has been determined numcri- 
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cally []J for g -> 0: 



n(0)A2 



rirX 



log 



380ft2 
mg 



(17) 



We can introduce (by hand) the concept of a critical den- 
sity into mean-field theory by selecting among the solu- 
tions t{jl) of eq. p4|) the one satisfying eg. pT]) . For 
the interaction parameters used in Ref [9|, mg/h^ = 



0.13, wc find a mean- field critical temperature t'^-^ 



Tkt /^bec = 0-69. This value is in excellent agreement 
with the Monte Carlo data. The inset of fig. [3] shows 
the variation of this mean-field critical temperature as a 
function of Cjz (for rag/h? ~ 0.13). 

The calculation of the mean-field critical temperature 
simplifies further in the strictly two-dimensional Bose 
gas, because the chemical potential in eq. (|13|) is then an 
explicit function of the critical density, and can be en- 
tered into eq. ([12]). With /ic = ln{l - exp (-7icA^)} ~ 
—mg/{?)8Qh^), and by again transforming the equation 
for N vs. density into a relation between critical temper- 
atures [^, [l^ , we obtain 



^2d 
^KT 



Ti2d 
rp2d 

-'bec 



Smg 
1^ 



In 



380^2 
mg 



-1/2 



(18) 



where we have neglected small corrections of order 
/2c log I /ic I- The strictly two-dimensional limit t^rj, = 
0.745 for mg/h'^ = 0.13 agrees with the data shown in 
the inset of fig. [3] in the large uJz limit. 

In conclusion, we have considered in this letter the 
semiclassical description of the quasi two-dimensional 
trapped Bose gas. We have compared this description 
with Quantum Monte Carlo data and have shown that 
the density profiles are accurately reproduced in the nor- 
mal phase down to the Kosterlitz-Thouless temperature. 
The thermal occupation of excited states in the out-of- 
plane direction quantitatively explains the large devia- 
tions of the density profiles from strictly two-dimensional 
mean-field theory, which was recently noticed in the ex- 
periment Originally, it was speculated that the 
emerging almost- Gaussian density profiles could be at- 
tributed to effects beyond mean-field 0, S|. However, 
even though the thermal occupation of the excited states 
in the tightly confined direction is clearly noticeable for 
the experimental parameters, the transition itself is still 
of the Kosterlitz-Thouless type as revealed by the exper- 
imental coherence patterns [J] and confirmed by numer- 



ical calculations of the algebraic decay of the condensate 
density with system size [9|. 
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